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The electronic properties of the polar interface 
between insulating oxides is a subject of great 
current interest [1 3j. An exciting new develop- 
ment is the observation of robust magnetism [4r 
[8] at the interface of two non-magnetic materi- 
als LaAlOa (LAO) and SrTiOa (STO). Here we 
present a microscopic theory for the formation 
and interaction of local moments, which depends 
on essential features of the LAO /STO interface. 
We show that correlation-induced moments arise 
due to interfacial splitting of orbital degeneracy. 
We find that gate-tunable Rashba spin-orbit cou- 
pling at the interface influences the exchange in- 
teraction mediated by conduction electrons. We 
predict that the zero-field ground state is a long- 
wavelength spiral and show that its evolution in 
an external field accounts semi-quantitatively for 
torque magnetometry data [5j. Our theory de- 
scribes qualitative aspects of the scanning SQUID 
measurements [6j and makes several testable pre- 
dictions for future experiments. 

Recent experiments on the LAO/STO interface have 
seen tantahzing magnetic signals [4l-(8], often persisting 
up to high temperatures ^ 100 K. A large magneti- 
zation of 0.3 — 0.4/iB per interface Ti was observed by 
torque measurements [S] in an external field. In con- 
trast, scanning SQUID experiments |5] found an inho- 
mogenous state with a dense set of local moments, but 
no net magnetization in most of the sample except for iso- 
lated micron-scale ferromagnetic patches. Our goal is to 
reconcile these seemingly contradictory observations and 
to gain insight into the itinerant versus local moment na- 
ture of the magnetism, the exchange mechanism, and the 
ordered state. 

LAO and STO are both band insulators, but the Ti02 
layers at the interface are n-doped when LAO is termi- 
nated by a LaO+ layer. The "polar catastrophe" jlj aris- 
ing from a stack of charged LaO+ and AlOj' layers grown 
on STO is averted by the the transfer of 0.5 electrons 
per interface Ti. In addition, oxygen vacancies are also 
known to provide additional electrons at the interface ^ . 

What is the fate of these electrons at the interface? 
Transport data suggests that only a small fraction of the 
electrons (5-10% of the 0.5 e/Ti) are mobile, as estimated 
from the Hall effect and quantum oscillations JUJ. 
Interestingly, magneto-transport studies show a large, 
and gate-tunable, Rashba spin-orbit coupling (SOC) for 



the conduction electrons, arising from broken inversion 
at the interface [TI]. In fact, most of the electrons (com- 
parable to 0.5 e/Ti) seem to behave like local moments 
in the magnetic measurements [3 [6] discussed above. 

In this paper we propose a microscopic model of elec- 
trons in the Ti t2g states at the LAO/STO interface and 
show that it leads to the following results. (1) 5=1/2 lo- 
cal moments form in the top Ti02 layer due to Coulomb 
correlations, with interfacial splitting of t2g degeneracy 
playing a critical role. (2) Conduction electrons medi- 
ate ferromagnetic interactions between the moments via 
double-exchange. (3) Rashba SOC for the conduction 
electrons leads to a Dzyaloshinski-Moriya (DM) interac- 
tion and a "compass" anisotropy term with a definite 
ratio of their strengths. (4) The zero field ground state 
is a long wavelength spiral with a SOC-dependent pitch. 
(5) The spiral transforms into a ferromagnetic state in 
an external field H. 

We provide a semi-quantitative understanding of the 
torque magnetometry results for the magnetization 
M{H). We also reconcile some of the key differences be- 
tween the torque and SQUID measurements and point to 
a novel magnetoelastic coupling effect that can be crucial 
for determining the magnetic ground state of this kind of 
polar interface. Our model naturally explains the coex- 
istence of magnetism and superconductivity at very low 
temperatures [51 IHl IB]- Finally, we make a number of 
specific predictions that can be tested experimentally. 

Symmetry-based considerations: Many of our 
conclusions regarding magnetism can be understood 
qualitatively based on symmetry, provided we have lo- 
cal moments Sr on a square lattice in the x-y plane at 
the interface. The detailed microscopic analysis given 
below shows how such a square lattice is formed and 
also gives quantitative insights into the parameter de- 
pendence of exchange couplings. Symmetry dictates the 
form of the interactions, the first of which is an isotropic 
Heisenberg exchange — JX^r ■ Sr+t/i, where (1 — x,y 
and h the lattice spacing. The sign of J is not deter- 
mined by symmetry, but our microscopic analysis leads 
to a ferromagnetic J > 0. Inversion symmetry breaking 
at the interface implies that ±z are not equivalent and 
we can write down two SOC terms [T2]. The DM term 
is D ^ (i^-(SrXSr+f,/i) with = zxfi, and the al- 
lowed compass anisotropy terms are — J2r\--^' i^r ^r+bx + 

One can see, quite generally, that the ground state of 
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such a model can be a long wavelength spiral that looks 
locally ferromagnetic (FM) , thus minimizing the J-term, 
but whose pitch is determined by the small DM and com- 
pass terms. Our microscopic results for D, A, A' unequiv- 
ocally predict such a spiral ground state for H = 0. In 
an external magnetic field (> IT), however, we find a 
uniform FM state. These observations permit us to un- 
derstand the torque data in a field semi-quantitatively, 
and to see why over most of the sample the H = scan- 
ning SQUID measurements find no net moment [Sj. We 
will discuss below detailed comparison with experiments 
and new predictions. 

Electronic structure: Electrons at the interface 
are in Ti t2g states. Both density functional theory 
(DFT) [13] and spectroscopic measurements [TiHTB] . 
show that, in the top Ti02 layer [labeled "1" in Fig.[T|a)], 
dxy states have lower energy than dxz , dyz . In addition 
to z-confinement raising the xz,yz levels, the mismatch 
of in-plane LAO and STO lattice parameters leads to an 
out-of-plane distortion that lowers dxy Further, the dxy 
orbitals delocalize primarily in the x-y plane, and simi- 
larly for dxz and dyz- DFT [T71 HH] and ARPES experi- 
ments |19| find the in-plane hopping between dxy orbitals 
t~0.3 eV, while the out-of-plane t' ~t/30. The resulting 
XZ (yz) bands are quasi one-dimensional (ID), dispersing 
primarily along x (y), and confined along z. 

Local moments: To form moments, we need to local- 
ize charge, either through disorder or correlations. Cou- 
pling to classical phonons is often argued to lead to lo- 
calization, but quantum effects lead only to a bandwidth 
reduction; localization needs disorder or correlations even 
in this case. Experiments |20) show that even the quasi- 
ID xz,yz bands, which should be the most sensitive to 
disorder, continue to contribute to transport. This rules 
out strong Anderson localization. 

The key to understanding the effect of correlations is 
the splitting of the t2g degeneracy at the interface de- 
scribed above. The 0.5e/Ti give rise to a quarter-filled 
xy band in the top Ti02 layer. We find that a modest on- 
site Hubbard for U and next-neighbor Coulomb V then 
lead to a charge-ordered insulator (COI), see Fig. [ijb). 
We obtain a simple analytical result for the phase bound- 
ary between a metal and COI insulator using a slave-rotor 
approach [21] (see Methods) . We also show that coupling 
to the breathing mode phonon further stabilizes the COI. 
DFT calculations need a rather large C/ ~ 8 eV to sta- 
bilize a COI [22]. We see from Fig. 1 that, for realistic 
values [23] of [/ = 4 eV {U/t ~ 13) and F ~ 0.5 - 1 cV 
{V/t ~ 1.5 — 3), we are deep in a checkerboard COI state: 
one sublattice occupied and the other empty. 

Exchange: We have now established a COI with 
S — 1/2 local moments. As argued in the introduction, 
symmetry then dictates the form of the effective mag- 
netic Hamiltonian HeS- While we can make progress at 
a phenomenological level, it is useful to ask what further 
insights a microscopic analysis leads to. 
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FIG. 1: Electronic orbitals and charge localization: (a) 

Schematic of Ti t^g energy levels. At the interface (layer 
1), the degeneracy is lifted with dxy lower than the dyz and 
dxz orbitals. The t2g energies in layers 2 and 3 increases 
due to the confining potential Konf- (b) Phase diagram of 
the single band extended Hubbard model at quarter-filing. 
A charge-ordered insulator (COI) is obtained for moderate 
values of correlations, on-site Hubbard U and next-neighbor 
Coulomb V, above the solid red phase boundary. Coupling to 
the breathing mode phonon further stabilizes the COI. The 
dashed black line corresponds to an energy gain Uph = 0.25t 
due to lattice distortion (see Methods). 



Consider then the problem of local moments in dxy 
orbitals on a checkerboard lattice [Fig. [2]^ a)] interacting 
via Hund's coupling Jh with a small density ric of con- 
duction electrons in xz, yz bands (possibly coming from 
doping by oxygen vacancies). Given Jh — 1 eV, it is 
much more reasonable to work in the non-perturbative 
double-exchange limit Jjj 2> i than in the perturbative 
RKKY hmit Jh < t. We treat the S" = 1/2 moments 
classically. This is justihed a posteriori because of the 
small quantum fluctuations in the magnetic states - FM 
or long-wavelength spiral - that we analyze. 

We find Hcff = Er i'Uxir) + Hyir)] + W with 



T~Lx{y) — — JSr-Sr-(-2£ ' 



A QV QV 



DyiS,xS,+2£) (1) 



where r denotes positions of occupied sites in the COI; see 
Methods. T-Lyir) is obtained by interchanging x ^ y and 
replacing -Din •H^(r). The term W = - J' J2r 
{Sr+x+y + Sr+s-y) couplcs the nearest neighbour spins 
along the diagonal direction [Fig. [2ja)]. The spins have 
normalization = 1. The form of Hx and Hy, as well 
as J ^ J', arises from the two-sublattice structure of the 
COI and the quasi-lD nature of xz, yz bands. We find 
J ~ rict/A and J' ~ Uct' /A for low carrier density Uc (the 
number of electrons per Ti in each band). 

The novel aspect of this double exchange model is 
the Rashba SOC 2Xaz • (k x cr) for the conduction elec- 
trons, where k is their momentum and cr are Pauli ma- 
trices. The SOC strength A is determined by the elec- 
tric field at the interface and transport experiments jllj 
find a gate voltage tunable A ~ 1.3—6.3 meV, so that 
X/t ~ 0.004-0.022. SOC then gives rise to the DM and 
compass terms in "Heff via the double exchange mecha- 
nism. For X/t ^ 1, we find D ~ ricA, A ~ 2ncA^/t 
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FIG. 2: Spiral ground state and its evolution with field: 

(a) Local moments in the COI with FM exchange interac- 
tions J and J'. The anisotropy and DM terms couple the 
same spins as J. (b) The zero-field spiral ground state, whose 
evolution with H = H{sm6HX + cos 6hz) with 6h = 15° is 
shown next, (c) The //-dependence of the spiral wavevector 
Q normalized by its i/ = value Qo for SOC X/t — 7r/25. 
(d) In-plane (M^) and (e) out-of- plane (Af^) components of 
the magnetization obtained from T = variational calcula- 
tions and low temperature (T = 0.25J) Monte Carlo (MC) 
simulations. The //-axis is scaled by //peak at which Mx has 
a maximum. Panels (f), (g), and (h) show the MC spin struc- 
ture factor /(q). (f) At // = the peaks of the MC /(q) lie on 
a white circle, the locus of wavevectors Qo corresponding to 
degenerate variational ground states, (g) For < H < //peak, 
the /(q) shows a net FM moment in addition to spiral peaks, 
(h) A FM state is obtained for H > //peak- 



and A' = 0. Two points are worth noting. The ratio 
AJ/D^ = 1/2 is independent of all microscopic parame- 
ters. Remarkably, our double-exchange results D/J^X/t 
and A/ J ^ [X/tY' are in the end very similar to Moriya's 
results [12] for superexchange. 

Spiral ground state: We have examined the ground 
state properties of "Hog in an external field H using a 
variety of analytical and numerical techniques. Let us 
first discuss H^O. Using a variational approach we find 
that for AJ/D^ < 1 the zero-field ground state is a spiral 



Sr — sin(Qo • r)Qo + cos(Qo • r)z. 



(2) 



with Qo = {2 X / ta) (cos if x + sin (py) where a is the Ti-Ti 
distance; see Methods. For AJ/D^ > 1, the ground state 
is an in-plane FM with Sr = cosi^ai; -I- sin 1^9 y. For all 
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FIG. 3: Torque magnetometry: (a) Magetization Mr = 
t/H (black symbols) from torque experiments i5, at T = 300 
mK with 6h = 15° shown in inset; (see Supplementary In- 
formation for background subtraction). Theoretical Mr{H) 
curve (red dashed line) computed with microscopic value 
AJ/D^ = 1/2 with y4~0.3 T chosen to match //peak- Better 
agreement (blue line) with experiment is obtained by the phe- 
nomenological choice AJ/D^ = 0.8 with A ~ 0.7 T; (see text). 
Inset shows schematic of two regimes in torque: an in-plane 
net magnetization M for small H evolving into a out-of-plane 
M for large H. (b) //peak scales with (X/t)^ for fixed ric, as 
predicted by the microscopic model, (c) Mt{H) curves, for 
different values of SOC that can be tuned by gating, collapse 
onto a single curve when H is scaled by //peak- 



AJ/D^ , the ground state is infinitely degenerate (with 
arbitrary < </? < 27r), even though HcS has Z4 symmetry. 
This peculiar degeneracy of compass anisotropy [23] is 
expected to be broken by "order-by-disorder" . 

The microscopically derived value AJ/D^ — 1/2 im- 
plies a zero-field spiral variational ground state. We have 
verified the stability of the spiral over FM at i/ = 
using two independent, unbiased numerical calculations 
(see Supplementary Information): finite-temperature 
MC and T ~ conjugate gradient energy minimization. 
In Fig. we show an example of spiral ordering in real 
space with a wavelength of a few lattice spacings. For re- 
alistic values of A/i ~ 0.02, our variational results give a 
pitch of (27r/Qo) - 600 A. 

In Figs. [2f^c)-(h), we show the evolution of the spiral 
state for X/t = tt/25 as a function of H applied in the 
xz-plane at an angle On = 15° with z-axis, (the geometry 
used in the torque magnetometry experiments [S]). In the 
variational calculation [Fig. |2jc)] the spiral state (Q 7^ 
0) develops into a fully magnetized state above a field 
Hpcak, which corresponds to a maximum in the in-plane 
component of the magnetization M = (Sr). We plot 
in Figs. [2]jd,e) and Mz as functions of H, where we 
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also see excellent agreement between the variational and 
the unbiased MC results. Here we used larger SOC that 
led to smaller pitch spirals for the MC simulations. Below 
we use the variational approach to study long wavelength 
spirals for realistic SOC A. 

The evolution from the spiral state to the FM as a 
function of field is best seen in the MC spin structure 
factor /q oc |(Sq)p shown in Figs. [2]^f)-(h), where Sq is 
the Fourier transform of Sr. The coplanar spiral a,t H = 
shows two peaks in /q at q = ±Qoy [Fig. [2|f)], which 
lie on the white circle, the locus of Qo corresponding to 
the degenerate ground states in eq. At intermediate 
fields, a uniform component develops along with the non- 
zero q peaks [Figj2jg)], while for H > iJpcak one obtains 
a FM state [Fig.[2pi)]. 

Torque and magnetization: The torque r = MxH 
is non-zero only in the presence of spin-space anisotropy, 
otherwise M would simply align with H. In Fig. |3][a), 
we compare our variational results with the magnetiza- 
tion Mr = t/H vs. H derived from torque magnetome- 
try [5] , where is the component of M perpendicular to 
H. (The experimental data have been shown after back- 
ground subtraction; see Supplementary Information.) 

We can understand Mr{H) as follows. In our theory, 
an external magnetic field with an in-plane component 
{9h 7^ 0) first induces a uniform magnetization in the 
spiral state [Fig. [2 . This leads to an increasing for 
H < -ffpoak in Fig. Sla), where iJpcak, at which Mt is a 
maximum, depends on 6h- For H > -ffpeak, the ground 
state is a FM [Fig. [2], in which the DM term is irrelevant, 
but the compass term gives rise to a large easy-plane 
anistropy A ~ 0.15 - 0.3 T (with A ~ 0.016 - 0.022 and 
Tic — 0.05) that tries to keep M in the plane. For large 
enough H the field dominates over the anisotropy, pulls 
M out of the plane and Mr decreases. These two regimes 
are shown schematically in Fig.|3][a) inset. 

In Fig. [3](a), we show two theoretical M{H) curves 
which differ in their choice of AJ /D^ . The red curve cor- 
responds to the Rashba value AJjD^ — 1/2, with A — 
2nc}?lt = 0.3 tuned to match i/poak (with Xjt = 0.022 
and Tic = 0.05). The high field drop in Mr{H) is much 
too rapid compared with the data. In the blue curve, we 
choose AJ/D'^ = 0.8 with X/t = 0.028 which leads to 
better agreement with experiment. Such a phenomeno- 
logical choice of AJ/D^ amounts to changing the planar 
anisotropy while keeping the DM interaction fixed. In our 
microscopic theory, both A and D were determined by 
Rashba SOC. In addition, several other effects also con- 
tribute to A, including dipolar interactions and atomic 
SOC (see Supplementary Information). 

The peak position Jfpcak oc ricX^ /t in our microscopic 
model and Mr{H) for different A collapse onto a single 
curve when plotted versus H / H^^ak', see Figs. [3]jb),(c). In 
practice, other sources of anisotropy might modify this 
scaling. Nevertheless, both ric and A are gate tunable 
and we expect iJpcak to vary substantially with bias. 



Ferromagnetic patches: We now discuss how we 
can reconcile the torque data [5] discussed above with the 
scanning SQUID results of ref. [6]. A detailed modeling of 
inhomogeneity seen in ref. [S] is beyond the scope of this 
paper, nevertheless, we can offer a plausible picture based 
on our theory. At iJ = most of the sample has, in fact, 
coplanar spiral order and hence no net magnetization. 
However, the energy gain of the spiral over FM order, 
Ae ~ [D"^ / J - A) /2^n^\^ /t'^Q.l-Q.2 K, is quite small; 
see Methods. Therefore small terms ignored up to now 
might well upset the balance in favor of FM. In particular, 
micron size patches could be possible if the coupling to 
strain fields was involved. 

We thus look at the effect of magnetoelastic coupling 
associated with polar distortion, a novel aspect of these 
interfaces. The Rashba SOC \{u) is determined by local 
electric fields, which are affected by small displacements 
uz of oxygen ions that bridge the Ti-O-Ti bonds [25]. 
Both A and D terms in HeS try to maximize energy 
gain by increasing \{u) by a displacement u < Q with 
A(m) = A— Aiw where Ai > 0; see Methods. But this costs 
an electrostatic energy —qEu, where E = —Ez is the 
electric field acting on the charge q of the oxygen ion, in 
addition to an elastic cost Kv? /2 with elastic constant 
K . Minimizing the total energy with respect to u, we find 
that the balance between spiral and FM can be reversed; 
see Methods. We hypothesize that strain influences the 
microscopic parameters in such a way that it stabilizes 
FM patches, which can then point in any direction in the 
plane, consistent with compass anisotropy. 

Discussion: We comment briefly on how our work, 
which builds on insights from electronic structure cal- 
culations [13l [T7J [iHl |22] , differs from other theoretical 
approaches. It is hard to see how one can get the large 
exchange coupling or net moment seen in LAO/STO from 
itinerant models [17] , which may well be relevant for 
other interfaces like GdTiOa/SrTiOa. Our model differs 
from ref. [28] in several respects. We are in the non- 
perturbative double-exchange limit, as noted above, and 
not the RKKY regime that leads to much smaller ex- 
change. More importantly, in our model, xz,yz carriers 
mediate exchange. The xy electrons in lower layers have 
essentially no interaction with moments in the xy orbitals 
of the top Ti02 layer. This results from the small overlap 
t' between xy orbitals along z as well as the level mis- 
match due to confinement. Conventional superconduc- 
tivity arises in xy states in lower layers and is decoupled 
from the local moments. 

Conclusions: Our theory has several testable predic- 
tions. Local moment formation is associated with (tt, tt) 
charge order at the interface. The zero-field state is pre- 
dicted to be a coplanar spiral state with several strik- 
ing properties, including a wavevector that scales linearly 
with Rashba SOC and hence is gate tunable. The evo- 
lution from a spiral to a FM state in an external field 
has characteristic signatures in the spin structure factor. 
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The easy-plane anisotropy strength is also dominated by 
SOC and exhibits substantial gate voltage dependence 
that can be tested in torque magnetometry experiments. 
The FM exchange J = n^t/A ~ 40 K (for = 0.05 
and t = 0.3 eV) should also be tunable by changing ric, 
the density of carriers. A theoretical analysis of the fi- 
nite temperature properties of our model is left for future 
investigations. 
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the support of NSF-DMR-1006532 (O.E., M.R.) and 
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Methods 

Effective Spin Hamiltonian: Broken inversion 
along z and in-plane square-lattice symmetry constrains 
the terms arising from SOC. In (ip-(SrXSr+6/i) the 

only allowed DM vector is = zxji. Other choices like 



dfj, = z break 7r/2 rotation and dx=±x,dy = ±y break in- 
plane reflection. Similarly, only the compass terms shown 
in the text are permitted. Terms that involve S'^S'^ or 
5^5*^ or S^S^ are not allowed by symmetry. 

Charge ordering Mott transition: We consider the 
quarter-filled, extended Hubbard model on a 2D square 
lattice with onsite U and nearest-neighbor Coulomb V. 
Using slave-rotor |2T] mean field theory (MET) for U 
and Hartree-Eock for V, we obtain an analytical result 
for the transition from metal to charge ordered insulator 
(COI). The checkerboard COI is stable for U > 8{t) when 
U < W, and for {U - 2V)V/U > (t) when U > W. 
Here (t) ~ 0.66t is the kinetic energy of occupied states. 
Slave rotor MET has been found to be in excellent semi- 
quantitative agreement with dynamical mean field theory 
(DMET) for several problems. In fact, our 2D square 
lattice results are quite similar to DMET of the Wigner- 
Mott transition on a Bethe-lattice [29] . 

The COI is stabilized by coupling to the breathing 
mode phonon, which has the same (tt, tt) wavevector as 
the charge order; see Supplementary Information. The 
only change above is that V ^ V = V + Eph, where Eph 
is the energy gained though lattice distortion. Typical 
values of Eph ~ 0.05 - O.lOeV [30]. 

Double Exchange: The local moments Sr on the 
sites of a 2D checkerboard lattice are described by their 
orientation {9r,^Pr)- In the large Jh limit, a conduc- 
tion electron on the "local-moment sublattice" (a) has its 
spin parallel to Sr. Thus we write the spin- full fermion 
operators draa-, with orbital index a = {xz),{yz), in 
terms of spinless fermions 0^0,, via at^„ — >■ F^^aJ^. Here 



Fr^==cos(6ir/2) and = sin(6'r/2)e 

Both spin projections a =ti4- are allowed on "empty 
sublattice" [h] sites, for which we use a common (global) 



quantization axis. The kinetic energy is then given by 

^A' ^ ^ ^ ^a/i -^rcr'^ra (^r+/i,aa- ^" ^r— /i,aa-) ^" h.C. (3) 



where i(a;2),s — ^(yz),^ — * and i(a;2),t; — — — t/ZQ. 
The Rashba terms tr^ sin(fcj,a) ^ ~kya{cr x z)y and 
sv[\{kxa) ^ kxa{<T x z)^ lead to the SOC Hamiltonian 

We can rewrite this in terms of the spinless a's as 



Hfl = ^ ^aaFrsa\jyh^+l,^o,a - br-il,acr) + h.C. (5) 
raa 



where a — —a and we define X(^xz)'\- = ~X{xz)i = and 

To obtain the parameters J, J', D and A of Hcs of 
eq. ([T]) starting from the microscopic "Hde = 'Hk+'Hr, we 
match the the energies of several low-lying configurations 
computed for both Hue and HcS- See Supplementary 
Information for details. 

Variational calculation: We study the ground state 
properties of HcS using a variational ansatz Sr = M -|- 
R(Q) cos(Q • r) + I(Q) sin(Q • r) with Q 7^ 0. The 
normalization S^ = 1 is satisfied via the constraints 
M • R = M • I = R • I = 0; R2 I2 and M2 + R2 = 1. 
We numerically minimize the energy per spin e to obtain 
optimal values for M, R, I and Q as a function of H, 
and calculate the torque. 

At zero field, we find a spiral ground state of eq. ^ 
with Q = Qo and M = for AJ / D'^ < 1 and an in- 
plane EM for AJ/D^ > 1. We can see this analyti- 
cally using the ansatz Sr — Rcos(Q • r) + Isin(Q • r). 
We write R = sin 9 sin ipx — sin 9 cos ipy + cos 9z and 
I — cosfx + sin (py with 9 e [0,7r] and (p £ [0, 27r). We 
first minimize with respect to Q (for Qa <C 1), find the 
optimal Q = {D /2 Ja) cos 9{cos ipx + sincpy) and the en- 
ergy e(Q) ~ -2{J+J')-A+{A~D^/J){cos^9)/2, which 
is independent of (p. For AJ/D^ < 1, the optimal choice 
of cos 9 = 1 leads to the spiral of eq. ^ with Q = Qo 
and an energy gain Ae = {D^/J — A)/2 relative to the 
EM state. For our microscopic model Ae = ricX^/t. 

The effect of a polar distortion uz of oxygen ions on 
the SOC can be modeled by replacing A — > A{u) = 
2ncA^(u)/t and I? — > D{u) = ncX{u) in e(Q) above, in 
addition to the energy cost [Kv? /2 — qEu). Here K is 
the elastic constant for out of plane distortion at the in- 
terface and E = —Ez is the electric field acting on oxygen 
ions at the interface. A(m) = A— Aiu with Ai > since the 
local electric field, and hence SOC, decreases if negatively 
charged oxygen ion moves up (u>0). We minimize the 
total energy with respect to u to find that energy gain 
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of the spiral is reduced: Ae ~ (ncA^/i)[l — 2qE\i/ KX]. 
With suitable choice of parameters, we can make Ae<0 
and stabilize the FM state over the spiral. 



[I] Ohtomo, A. & Hwang, H. Y. A high-mobility electron gas 
at the LaAlOs/SrTiOs heterointerface. Nature 427, 423- 
426 (2004). 

[2] Thiol, S., Hammcrl, G., Schmehl, A. Schneider, C. W. & 
Mannhart, J. Tunable quasi two-dimensional electron 
gases in oxide heterostructures. Science 313, 1942-1945 

(2006) . 

[3] Zubko, P., Gariglio, S., Gabay, M., Ghosez P. & Triscone, 
J.-M. Interface physics in complex oxide heterostructures. 
Annu. Rev. Condens. Matter Phys. 2 141-165 (2011). 

[4] Brinkman, A et al. Magnetic effects at the interface be- 
tween non-magnetic oxides. Nature Materials 6, 493-496 

(2007) . 

[5] Li, L., Richter, C, Mannhart, J. & Ashoori, R. C. Coexis- 
tence of Magnetic Order and Two-dimensional Supercon- 
ductivity at LaAlOs/SrTiOs Interfaces. Nature Phys. 7, 
762-766 (2011). 

[6] Bert, J. et al. Coexistence of Magnetic Order and Two- 
dimensional Superconductivity at LaAlOa/SrTiOs Inter- 
faces. Nature Phys. 7, 767-771 (2011). 

[7] Ariando et al. Electronic phase separation at the 
LaAlOs/SrTiOa interface. Nature Comm. 2, 188-194 
(2011). 

[8] Dikin, D. A., Mehta, M., Bark, C. W., Folkman, CM., 
Eom, C. B. & Chandrasekhar, V. Coexistence of super- 
conductivity and ferromagnetism in two dimensions. Phys. 
Rev. Lett. 107, 56802-56805 (2011). 

[9] Shalom, M. B., Ron, A., Palcvski, A. & Dagan, Y. 
Shubnikov-De Haas oscillations in SrTiOa/LaAlOa inter- 
face. Phys. Rev. Lett. 105, 206401-206404 (2010). 

[10] Lerer, S., Shalom, M. B., Deutscher, G. & Dagan, 
Y. Low-temperature dependence of the thermomagnetic 
transport properties of the SrTiOa/LaAlOs interface. 
Phys. Rev. B 84, 075423-075428 (2011). 

[II] Caviglia, A. D., Gabay, M., Gariglio, S., Rcyren, N., 
Cancellieri, C. & Triscone, J.-M. Tunable Rashba spin- 
orbit interaction at oxide interfaces. Phys. Rev. Lett. 104, 
126803-126806 (2010). 

[12] Moriya, T. Anisotropic superexchange interaction and 
weak ferromagnetism. Phys. Rev. 120, 91-98 (1960). 

[13] Pentcheva, R. & Pickett, W. E., Charge localization or 
itineracy at LaAlOs/SrTiOs interfaces: Hole polarons, 
oxygen vacancies, and mobile electrons. Phys. Rev. B 74, 
035112-035118 (2006). 

[14] Salluzzo, M. et al. Orbital reconstruction and the two- 
dimensional electron gas at the LaAlOa/SrTiOa interface. 
Phys. Rev. Lett. 102, 166804-166807 (2009). 

[15] Sing, M. et al. Profiling the interface electron gas of 
LaAlOa/SrTiOs heterostructures with hard x-ray pho- 
toelectron spectroscopy. Phys. Rev. Lett. 102, 176805- 
176808 (2009). 

[16] Berner, G. et al. LaAlOa/SrTiOa oxide heterostructures 
studied by resonant inelastic x-ray scattering. Phys. Rev. 
B82, 241405-241403(R) (2010). 

[17] Popovic, Z. S., Satpathy, S. & Martin, R. M. Origin of the 
two-dimensional electron gas carrier density at the LaAlOs 



on SrTiOs interface. Phys. Rev. Lett. 101, 256801-256804 

(2008). 

[18] Hirayama, M., Miyake, T. & Imada, M. Ah initio Low- 
energy model of transition-metal-oxide heterostructure 
LaAlOs/SrTiOa. J. Phys. Sac. Jpn. 81, 084708-084716 
(2012). 

[19] Santander-Syro, A. F. et al. Two-dimensional electron 
gas with universal subbands at the surface of SrTiOs. Na- 
ture 469, 189-194 (2011). 

[20] Fete, A., Gariglio, S., Cavigha, A. D., Triscone, J.-M. & 
Gabay, M. Rashba induced magnetoconductance oscilla- 
tions in the LaAlOs-SrTiOa heterostructure. Phys. Rev. B 
86, 201105-201108(R) (2012). 

[21] Florens, S. & Georges, A. Slave-rotor mean-field theories 
of strongly correlated systems and the Mott transition in 
finite dimensions. Phys. Rev. B 70, 035114-035128 (2004). 

[22] Pentcheva, R. & Pickett, W. E., Ionic relaxation con- 
tribution to the electronic reconstruction at the n-type 
LaAlOa/SrTiOs interface. Phys. Rev. B 78, 205106- 
205110 (2008). 

[23] Imada, M., Fujimori, A. & Tokura, Y. Metal-insulator 
transitions. Rev. Mod. Phys., 70, 1039-1263 (1998). 

[24] Nussinov, Z., Biskup, M., Chayes, L. & van den Brink, J. 
Orbital order in classical models of transition-metal com- 
pounds. European Phys. Lett. 67, 990996 (2004). 

[25] Khalsa, G., Lee, B. & MacDonald, A. H. Theory of t2g 
electron-gas Rashba interactions. arXiv 1301.2784 (2013). 

[26] Chen, G & Balcnts, L. Ferromagnetism in itinerant two- 
dimensional t2g systems. arXiv 1301.5681 (2013). 

[27] Fischer, M. H., Raghu, S. & Kim, E.-A. Spin-orbit cou- 
pling in LaAlOa/SrTiOa interfaces: magnetism and or- 
bital ordering. New J. Phys. 15, 023022-023032 (2013). 

[28] Michacli, K, Potter, A. C. & Lee P. A. Superconducting 
and ferromagnetic phases in SrTiOa/LaAlOa oxide inter- 
face structures: possibility of finite momentum pairing. 
Phys. Rev. Lett. 108, 117003-117007 (2012). 

[29] Camjayi, A., Haule, K., Dobrosavljevic & V., Kotliar, 
G. Coulomb correlations and the Wigner-Mott transition. 
Nature Phys. 4, 932-935 (2008). 

[30] Nanda, B. R. K. &: Satpathy, S. Electronic phases and 
phase separation in the Hubbard-Holstcin model of a polar 
interface. Phys. Rev. B 83, 195114-195120 (2011). 



7 



SUPPLEMENTARY INFORMATION 

1) Electron-phonon coupling in charge-ordered 
insulating state: Here we show how the charge-ordered 
insulator (COI) is further stabilized by coupling to the 
breathing mode phonon. The electron-phonon cou- 
pling can be written as Hep = {Kq/2) Y^iiQ^ai + Qli) - 
9 J^iiQaii^ai + Qbinbi)- Here a, b label the two sublattices 
of the COI in the i-th unit cell, Kq is the elastic constant, 
g the electron-phonon coupling strength for breathing 
model displacements Q's of Ti ions, and the 
electronic occupancies. In the COI state with one sub- 
lattice occupied Uai — 1 and the other empty — 0, we 
take the staggered displacement pattern Qai = —Qm = Q 
and minimize the energy with respect to Q. The energy 
gained though such lattice distortion is Eph = g'^ /AKi^), 
which has a typical value of ~ 0.05 — O.lOeV [T]. Incor- 
porating the effects of H^p into our slave rotor analysis, 
we find that the results are the same as those described 
in the Methods with the replacement V V = V + Eph- 
This enhancement in V is due to the cooperative distor- 
tion of the breathing mode as the charge order have the 
same (tt, tt) wavevector. 

2) Microscopic derivation of parameters for 

'Heff- We obtain the parameters J, J', D and A of Hoff by 
comparing the energies of several low-lying configurations 
with that obtained from the double exchange Hamilto- 
nian "Hde = "Hk + "Hr described in the Methods. Here 
we briefly sketch the procedure. 

Consider an iV x TV (with even N) lattice of Ti sites 
with local moments in a checkerboard arrangement 

(see Fig. 2(a) in the main text). Consider conduction 
electrons in the low-density limit nc ^ 1 so that only the 
states near the bottom of the xz and yz bands are filled. 

We first neglect t' <^ t and take strictly ID xz and 
yz bands. This leads to decoupled ID channels along 
X and y. For t' = 0, there is no coupling (J' = 0) 
between nearest-neighbor spins of the COI along diag- 
onal direction. The system then essentially splits into 
two inter-penetrating, decoupled subsystems. (We will 
consider the t' ^ coupling between these subsystems 
below.) Let us consider the energies of the following two 
configurations of local moments: (a) uniform FM state 
and (b) a twisted spin configuration. 

(a) An uniform FM with 9^ = and ip^ — ip. We 
diagonalize Hde to obtain the lowest-lying energy bands 
for xz and yz carriers, given by 

A 

CxzO^) — — 2i -I- t{kx — — sin 6 sin 0)^ — sin^ 6 sin^ (p, 

and eyz(k) obtained by replacing ky and sintp 

cos (fi. Depending on the spin configuration the system 
gains energy by shifting the band bottoms. This essen- 
tially leads to the easy-plane anisotropy as can be seen by 
computing the total energy per spin to leading order in 



ric and X/t, namely, euE{S, ^fi) — —Srict— {2nc}? /t) siv? 9. 
This can be compared with the energy of the same con- 
figuration in the spin- model, eofi(6', ip) ~ — 2J — Asit? 9 
to obtain the result A = 1x1^}? jt mentioned in the main 
text. (The angle-independent constant determines the 
absolute value of energy and is not a meaningful quantity 
to compare between microscopic and low-energy effective 
models.) 

(b) A twisted spin configuration along x direction with 
^'r+2x — 6'r = (J^* — 0, 9j.j^iy = 9r and (fir = 0. In this 
case the xz band gives rise to N identical ID channels 
at Uy = 1,..,A^. One can analytically diagonalize the 
xz part of the Hamiltonian, by transforming the &-site 
fermions into new rotated fermions /'s as 

br+x,xzt = COS (6'r/2) /r+i,a;^t ^ ^iu {9r/2) fr+x,xzl, (6a) 
br+x,xzl = e*'^'-[sin (9^/2) fr+x.,xz\ + COS {9j./2) fr+x.xzi], 

(6b) 

provided {N/2)59 < tt so that < 6'r < tt and i^r = 
Vr. The energy per spin from this part turns out to be 
exz{59) ~ -Anct-nc\59+{nJ,/^)59'^ . The yz band leads 
to N decoupled ID channels at Ux = 1,..,A^, but each 
with a separate FM configuration of spins with angles 
9 — nx59/2 and (/? = 0. Here we can use the result for 
eyz{\<i) obtained in (a) and calculate the energy contribu- 
tion per spin from yz band as €yz{59) ~ —Anct — 2nc\'^/t. 
The total energy for this configuration is eDE('50) — 
—^nct—2ncX?/t—ncX59 + {nct/%)59^ compared with that 
obtained from U^s, e^s{59) ~ -2J - A- D59 +{J /2)59'^ . 
This gives J = Uct/A and D = UcX. 

Finally, we switch on t' ^ and obtain J' (see 
Fig. 2(a) in the main text). Here we neglect the small 
SOC terms (again of the DM and compass form) that 
are suppressed by a factor {f /t) relative to D and A. 
We now consider a twisted spin configuration (c) where 

Or+x+i, - 6'r = S9, 9r+2x = Or+2y = Or and (fir = '-P ■ The 

fermion Hamiltonian can be diagonalized by rotating h- 
site fermions via equation ([6|, now for both xz and yz 
carriers. The energy in this case to leading order in t' jt is 
eDE(^6') ^ -%nc{t + t')^{jict' |^)b9'^ . The corresponding 
energy from spin model is eeff((50) — —23' + J' 59^ . This 
leads to a J' = Uct' /A < J. 

3) Monte Carlo and numerical minimization: 
We have performed Monte Carlo (MC) simulation [5] and 
conjugate gradient minimization [3] as unbiased checks 
on the results obtained from variational calculation for 
ground state properties of "Heff- In order to access spiral 
wavevector Qq = {2X/ta), we have taken commensurate 
values of A = n/p with p = 25, 30, 50 for 2xnpxnp spins 
{n = 1, 2) on the lattice of Fig. 2(a) (main text). We use 
standard Metropolis sampling [J' for MC simulation with 
unit length spins and perform 10^ MC steps per spin 
for equilibration, followed by 4 x 10^ MC steps/spin to 
calculate magnetization M = (Sr) and spin structure 
factor Jq cx |(Sq)p. The finite-temperature spin struc- 
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ture factor plotted in Figs. 2(f)-(h) is normalized with 
/ = X^BZ summed over the first Brillouin zone. 

4) Analysis of torque magnetometry data: In the 

main text [Fig. 3(a)] we compare our theoretical results 
with the background subtracted torque magnetometry 
data Here, we describe the background subtraction 
procedure used for the data and its rationale. We also 
discuss the scaling of the magnitude of the magnetization. 

Apart from M x H from the sample, the experimental 
torque r contains additional "background" contributions 
that are apparent at high fields _ff 5 T. First there 
is a paramagnetic piece of the torque Tqu ^ that 
originates from the STO substrate |3]. The raw data for 
the torque r and the measured paramagnetic signal Tqu 
from the STO substrate (extrapolated up to high fields) 
from ref. [4] are plotted in Supplementary Fig. |4]ja). 

We see from Supplementary Fig. |4|^b) that the sub- 
tracted torque (r — Tqu) is a linear function (rnn) of H at 
fields in excess of 10 T, which in fact has a zero (or negli- 
gible) intercept when extrapolated down to fl" = 0. Why 
is the high field behavior linear in H even beyond 15T? 
We consider two possibilities. One possible explanation 
could be an extremely high, and rather unprecedented, 
strength of easy plane anisotropy > 15 T. In this case, 
the magnetization M vector is not pulled out of the plane 
of the sample to completely align with the external field 
even at 15 T, for only then could one get a torque linear 
in H. This seems to us most unlikely; certainly none of 
the mechanisms we are aware of would account for such 
a large anisotropy. 

An alternative explanation (Lu Li, private communi- 
cation) is that the high field linear- signal arises from 
a magnetic field inhomogeneity that gives rise to force 
(M • V) H. Given the fact that the high-field part is 
found to vary from sample to sample [3], we think that 
this second explanation is much more viable. 

We have therefore subtracted this high-field linear 
background from the torque. In Fig. 3(a) (main text) 
we plot Mt{H) = (t — — Tnn)/H for the data of 
ref. P]. 

In comparing our theory with the experimental AIr{H) 
we have simply scaled our Mr(H-peak) to match the ex- 
perimental value in Fig. 3(a) (main text). However, we 
point out that there is quite reasonable agreement be- 
tween theory and experiment even if he had not scaled 
our theory and plotted the results in absolute units us- 
ing a gr-f actor of 2. We do not try and derive the exact 
scaling factor microscopically in view of the uncertain- 
ties in the actual ^-factors of the local moments (0.5 per 
Ti) and of the conduction electrons (2nc per Ti from the 
two bands) . Note that in the double-exchange model the 
conduction electron spins are locked to those of the local 
moments. 

5) Magnetic Anisotropy: In the text we focussed 
on anisotropy arising from Rashba SOC. Other mecha- 
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FIG. 4: Analysis of torque data: (a) Torque magnetometry 
data for LAO/STO interface from ref. [4]. Also shown is the 
quadratic background {jqu) measured for STO substrate, (b) 
Shows the high field Wneax-H contribution rn extrapolating 
through the origin. 



nisms include shape anisotropy (via dipolar interactions) 
and atomic SOC [5]. Although these are likely to be 
subdominant to the Rashba contribution, as argued be- 
low, nevertheless, the phenomenological fit in Fig. 3(a) 
(main text) suggests the importance of more easy-plane 
anisotropy than can be accounted for by Rashba SOC 
alone. 

FM domains of size ^ 1/im or larger would prefer to 
align in-plane in 2D or quasi-2D systems to minimize the 
long-range dipolar interaction. The energy contribution 
per unit volume due to such shape effect is [3 E] edip — 
(/io/2)m^ cos^ 6*, where m is the average magnetization 
per unit volume and 9 the angle it subtends with the 
normal to the plane of the interface. Phenomenologically, 
we can consider it as a term in the spin model of the 
form As X]r('^r whcrc the easy-plane anisotropy due to 
shape effect is Ag = ^q^jl^S"^ / Aa? , with S* = 1/2 and 2a? 
as the volume occupied by the individual local moments 
in the COL As a result. As ~ 0.015 T. Note that this 
crude estimate of shape anisotropy is much smaller than 
the compass anisotropy term A ~ 0.3 T that arises from 
Rashba SOC. 

Atomic SOC prefers the magnetization to point along 
crystalline axes, leading to four-fold in-plane anisotropy. 
This must be small since the FM patches, seen in scan- 
ning SQUID experiments [7] , point in random directions 
in the plane. 
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